Signatures of phase transitions in nuclei at finite 

excitation energies 

Y. Alhassid*, C. Ozen^ and H. Nakada** 

* Center for Theoretical Physics, Shane Physics Laboratory, Yale University, New Haven, CT 

06520, USA 

^Faculty of Engineering and Natural Sciences, Kadir Has University, Istanbul 34083, Turkey 
** Department of Physics, Graduate School of Science, Chiba University, Inage, Chiba 263-8522, 

Japan 

Abstract. The mean-field approximation predicts pairing and shape phase transitions in nuclei as a 
function of temperature or excitation energy. However, in the finite nucleus the singularities of these 
phase transitions are smoothed out by quantal and thermal fluctuations. An interesting question is 
whether signatures of these transitions survive despite the large fluctuations. The sheU model Monte 
Carlo (SMMC) approach enables us to calculate the statistical properties of nuclei beyond the mean- 
field approximation in model spaces that are many orders of magnitude larger than spaces that 
can be treated by conventional diagonalization methods. We have extended the SMMC method to 
heavy nuclei and used it to study the transition from vibrational (spherical) to rotational (deformed) 
nuclei in families of rare-earth isotopes. We have calculated collective enhancement factors of level 
densities as a function of excitation energy and found that the decay of the vibrational and rotational 
enhancements is well correlated with the pairing and shape phase transitions, respectively. 
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INTRODUCTION 

A fundamental relation describing the thermodynamics of a system at temperature T is 
the free energy 

F(r) = -rinz(r) , (i) 

where Z(7') is the partition function. For a system with a Hamiltonian H 

Z{T) = Txe'^l'^ . (2) 

In the presence of correlations, // is a many-body Hamiltonian, and it is difficult to cal- 
culate its partition function. A tractable approximation is the mean-field approximation, 
such as the Landau theory [1]. In the mean-field approximation, the free energy F[T, a] 
is evaluated as a function of certain order parameters cr and then minimized with re- 
spect to o to find the equilibrium configuration. Examples of order parameters are the 
quadrupole deformation in the Landau theory of nuclear shape transitions [2], and a 
complex pairing gap in the theory of superfluidity. 

However, in a finite-size system such as the atomic nucleus it is often important 
to take into account fluctuations of the order parameters. In general there are two 
types of fluctuations: (i) thermal or static fluctuations, and (ii) quantal fluctuations. The 



probability for a static fluctuation <j is given by P[g] oc e ^I^'*'']/^, and the partition 
function is approximated by 

Z(r) =y"D[a]e-^[^'*'l/^. (3) 

Quantal fluctuations become important at low temperatures and are included by consid- 
ering time-dependent fluctuations a = o{t) of the order parameters. 

A systematic approach to include fluctuations beyond the mean-field approximation 
is based on the Hubbard-Stratonovich (HS) transformation [3], in which the imaginary- 
time many-body propagator is described as a functional integral over one-body propa- 
gators of particles moving in external time-dependent auxiliary fields. The saddle-point 
approximation leads to a mean-field theory. In the static path approximation (SPA) [4], 
large-amplitude fluctuations of the static fields are integrated over, while in the SPA 
plus random phase approximation (RPA), small-amplitude time-dependent fluctuations 
are taken into account around each static fluctuation [5, 6]. In the auxiliary-field Monte 
Carlo (AFMC) approach, also known in the context of the nuclear shell model as the 
shell model Monte Carlo (SMMC) method [7, 8], all fluctuation are included by a 
Metropolis sampling of the auxiliary-field configurations. 

In a mean-field approximation the nucleus is predicted to undergo pairing phase 
transition and shape phase transition as a function of temperature or excitation energy. 
Shape transitions are also predicted as a function of proton or neutron number. The 
fluctuations smooth the singularities of the phase transitions and an interesting question 
is whether signatures of these transitions survive despite the large fluctuations. 

In the following we discuss the SMMC method and use it to study the emergence of 
collectivity in the configuration-interaction shell model approach [9, 10]. In particular, 
we focus on rotational collectivity associated with deformed nuclei and vibrational col- 
lectivity associate with spherical nuclei. We identify a thermal observable that can dis- 
tinguish between the different types of collectivity and use it to study the transition from 
spherical to deformed nuclei in families of rare-earth isotopes. We define a collective 
enhancement factor of level densities and demonstrate that the decay of vibrational and 
rotational collectivity is correlated with the pairing and shape transitions, respectively. 

THE SHELL MODEL MONTE CARLO (SMMC) APPROACH 

We first discuss the SMMC method in the framework of the canonical ensemble (at fixed 
temperature) and then describe how to translate its results to fixed excitation energy. 



Hubbard-Stratonovich transformation for the canonical ensemble and 

the SMMC method 

The Gibbs operator e^^^ at inverse temperature j8 = l/T is also the system's prop- 
agator in imaginary time /3 . The Trotter product of the propagator is obtained by di- 
viding the time interval [0,/3] into Nt time slices of length Aj8 each, and rewriting 



. Expressing the two-body interaction as a quadratic form in the 

one-body densities, the short-time propagator e~^^^ can be written as a integral over 

auxiliary fields a of one-body propagators described by a one-body Hamiltonian h(y 
that is linear in the one-body densities. A different set of fields C7(t„) is required for 
each time slice T„ = nAj8 {n— 1, . . . , A'f). The finite-time propagator is then written as a 
functional integral 

e-P" = J D[a]GaUa (4) 

over time-dependent auxiliary fields ct(t) with a Gaussian weight G^- Here Ua = 
g-Aphc{TNt) g-A/3/i(7(zi) jj^g many-particle propagator for a time-dependent one-body 
Hamiltonian haiT) that describes non-interacting particles moving in external fields 
cr(T). Relation (4) is known as the Hubbard-Stratonovich (HS) transformation [3]. 
The thermal expectation value of an observable O is calculated from 

^ Tr{Oe-P") ^ fD[a]W^^^{0)a 

where Wa — Gcy|Tri[/fj| is a positive-definite function, 4>(j = Tr:Ua/\T''rUa\ is the Monte 
Carlo sign function, and {0)(j = Tr{OU(y)/TrU(y is the thermal expectation value of 
the observable at a given configuration of the auxiliary fields o. In SMMC, we use a 
Metropolis algorithm to choose configurations that are distributed according to the 
weight function Wa- The thermal expectation value in (5) is then estimated from 

In the finite nucleus, it is necessary to calculate observables at fixed number of protons 
Z and neutrons A^. To that end, we use the canonical ensemble in which the traces in (5) 
and in the weight function Wa are calculated at fixed Z and A'^. The projection on a 
fixed number of particles can be expressed as a Fourier sum. For example, for A^^ single- 
particle orbitals, the canonical partition function Tr^i Uq for A particles is given by 

TtaUo = £ e-'-^^^Tr ie^^^e'^'-^Ua) , (7) 

^■s m=l 

where A is the particle-number operator, 0^ = Inm/Ns (m = 0, . . . , A^^) are quadrature 
points and /i is a real chemical potential [used to stabilize numerically the sum in (7)]. 
The trace on the right-hand side of (7) is a grand-canonical trace evaluated over the Fock 

space with all values of particle number and is given by det ^1 + e"^'"e^'^U(y^ , where U(j 

is the Ns X A^^ matrix representing Uo in the single-particle space. 

The SMMC method is particularly useful in the microscopic calculations of statistical 
properties of nuclei [11]. Many of the early applications were to medium-mass nuclei, 
and more recently the method was extended to heavy nuclei, overcoming a number of 
technical challenges [9]. 



M icrocanonical ensemble 



In the SMMC method we calculate observables at fixed temperature, i.e., in the 
canonical ensemble. However, nuclear properties are usually measured at a fixed energy 
E, necessitating the use of the microcanonical ensemble 5{E — H). The microcanonical 
ensemble is related to the canonical ensemble by an inverse Laplace transform 

/i'oo 
dpe^^e-^" . (8) 
-i'oo 

The inverse Laplace transform is numerically unstable but can be evaluated in the saddle- 
point approximation. 

An example is the level density p{E) at energy E whose relation to the canonical 
partition function is obtained by taking the trace of Eq. (8) 

/oo 
dPe^^ZiP) . (9) 
-CXD 

The average level density is obtained by evaluating (9) in the saddle -point approxima- 
tion. The saddle-point condition —d\nZ/dj5 — E determines /3 as a function of E and 
the level density is approximated by 

where S{E) and C are the canonical entropy and heat capacity, respectively. In SMMC, 
we calculate the energy £ = £'(/3) as a function of /3 and integrate the find the partition 

function, i.e., lnZ(/3) — InZ(O) — Jq dj5'E{P'). The entropy is then calculated from 
S{E) = lnZ(j8 ) -I- PE, while the heat capacity is found by taking a derivative C = dE /dT. 
The latter derivative has large statistical errors at low T . These errors can be reduced 
significantly by calculating the derivative inside the HS integral and taking into account 
correlated errors [12]. 



EMERGENCE OF COLLECTIVITY IN THE SHELL MODEL 

APPROACH 

Heavy nuclei are known to exhibit various types of collectivity, e.g., vibrational collec- 
tivity in spherical nuclei and rotational collectivity in well-deformed nuclei. Models such 
as the geometric Bohr Hamiltonian [13] and the interacting boson model [14] have been 
successful in describing these types of collectivity. However, a microscopic description 
in the framework of a truncated spherical shell model has been a major challenge. In 
particular, heavy nuclei can have large deformation that are more difficult to reproduce 
in a truncated spherical model space. 

The dimensionality of the many-particle model space that is required to describe 
such heavy deformed nuclei is many orders of magnitude beyond the capability of 
conventional diagonalization methods [9]. The SMMC method enables calculations 
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FIGURE 1. Left: (J^)?- as a function of temperature T for ^^^Dy. The solid circles with statistical errors 
are the SMMC results. The solid line is a fit to the rotational band model and the dash-dotted line is a fit 
to the vibrational band model [see Eq. (1 1)]. Right: total state density p of '^^Dy versus excitation energy 
Ey;. The SMMC state density (solid circles) is compared with experimental results: level counting at low 
energies (histograms), neutron resonance data (triangle) and data measured by the Oslo method [16] (open 
squares). Adapted from Ref. [9]. 



within such large model spaces. However, SMMC is a suitable method for calculating 
thermal and ground-state observables, but not for calculating detailed level schemes. 
This poses a new challenge in the framework of SMMC since the specific type of 
collectivity is usually identified through its signatures in spectroscopy. 

This difficulty can be overcome by calculating thermal observables that are sensitive 
to the specific type of collectivity. We have identified such a thermal observable in (J^)r, 
where J is the total angular momentum of the nucleus [9, 10]. At low temperatures T 
this observable is dominated by the ground-state band and for an even-even nucleus [15] 

, 30—^ -T vibrational band 

(J')r- {}-e-'--l'y (11) 



^2 



rotational band 



where £2+ is the excitation energy of the first 2+ level. Thus in the limit of sufficiently 
low temperatures a rotational nucleus is characterized by the linear temperature depen- 
dence of (J^)r, in contrast with a vibrational nucleus for which (J^)r exhibits a non- 
linear dependence on temperature. 



Rotational collectivity in heavy deformed nuclei 

We first discuss the example of a strongly deformed rare-earth nucleus, ^^^Dy [9]. The 
left panel of Fig. 1 shows (J^)r as a function of temperature T . The solid circles are the 
SMMC results and the vertical bars are the statistical errors. The solid and dash-dotted 
lines are, respectively, fits to the rotational- and vibrational-band formulas in Eq. (11). 
The agreement of the straight line fit with the SMMC results confirms the rotational 
character of ^^^Dy. Furthermore, from the slope of the straight line we can determine 
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FIGURE 2. (1^)7- as a function of temperature T for a family of even samarium isotopes '"'^^^'^^Sm. 
The circles with error bars are the SMMC results and the solid lines are obtained from Eq. (12) using a 
complete set of low-lying experimental energy levels and an experimentally determined BBF level density 
(except for '^^Sm for which only discrete levels are used). Adapted from Ref. [10]. 



the ground-State moment of inertia /g.g. = 35.5 ±3.3 MeV^ , in agreement with the 
experimental value of /gj,. = 3/£'2+ = 37.2 MeV^^ 

The total state density of ^^^Dy is shown in the right panel Fig. 1 as a function 
of excitation energy E^. The SMMC state density (solid circles) compares well with 
various experimental results. The dash-dotted line is the level density obtained in the 
Hartree-Fock-Bogoliubov (HFB) approximation. The SMMC density is enhanced when 
compared with the HFB density since it includes rotational bands that are built on top of 
intrinsic band heads. 



Transition from vibrational to rotational collectivity in heavy nuclei 

We have studied families of even-even samarium (^^^^^^"^Sm) and neodymium 
^144-I52jsj^-) isotopes [10]. These families are known to exhibit shape transitions from 
spherical to deformed nuclei as the number of neutrons increases from shell closure 
towards the mid-shell region. In the HFB approximation we find that ^^^Sm is spherical 
in its ground state (i.e., at T = 0). ^^'^Sm acquires a small deformation in its ground state 
while ^^^Sm is a well-deformed nucleus. ^^^Sm is known as an X(5) nucleus defining 
the phase transition between spherical and axially-deformed nuclei [17, 18]. Similarly 
i44Nd and ^^^Nd are spherical while ^"^^Nd acquires a small deformation that increases 
with the number of neutrons. 

Fig. 2 shows (J^)r versus temperature T for the samarium isotopes. We can clearly 
identify in the behavior of (J^)r a crossover from the vibrational nucleus ^"^^Sm to the 
rotational nucleus ^^^Sm. The solid hnes are extracted from experimental data using 

(J')r = ^ {^h[J^ + 1)(27, + l)e-^'/^ + £dE,PiE.) (J')£. e^^^^^ j , (12) 

where Z{T) = Y.f{2Ji + \)e^^'/^ + J^^dExp{Ex)e^^''/^ is the corresponding partition 
function. The energies Ei define a complete set of low-lying experimental energy levels 
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FIGURE 3. The enhancement factor K asa function of excitation energy Ex in even samarium isotopes 
148-I54jgjj^ The thin arrows indicate the excitation energies of the proton and neutron pairing transitions 
and the thick arrows correspond to the shape transitions. Adapted from Ref. [10]. 



with spin 7,, while p (E^) is a backshifted Bethe formula (BBF) level density extracted 
from experimental data (i.e, level counting at low energies and neutron resonance data at 
the neutron threshold energy). The BBF level density is used above a certain excitation 
energy Ej^/ below which a complete set of experimental discrete levels is known. We 
observe an overall good agreement between the SMMC and experimental results. 



Phase transitions and the collective enhancement factors 

Collective states lead to enhancement of level densities. Most available expressions 
for the collective enhancement factors are phenomenological [19]. The level density 
calculated in the finite-temperature HFB approximation only counts the intrinsic states. 
Therefore, we can define collective enhancement factor by ^ = PsMMC / Phfb , measuring 
the ratio between the SMMC total state density and the intrinsic HFB density. In Fig. 3 
we show the calculated K (on a logarithmic scale) as a function of excitation energy Ex 
for the family of samarium isotopes ^"^^"^^^Sm. 

In the HFB approximation we can identify pairing and shape phase transitions which 
are indicated by arrows in Fig. 3. There are two pairing transitions, one for protons and 
a second for neutrons. ^^^Sm is spherical in its ground state so that only the pairing 
transitions occur In such a spherical nucleus the collective enhancement originates in 
vibrational states, and we observe that K decays to ~ 1 above the pairing transitions. The 
heavier samarium isotopes shown in Fig. 3 are all deformed in their ground state and K 
has a minimum above the pairing transitions. These observations suggest that the decay 
of vibrational collectivity is correlated with the disappearance of pairing correlations. In 
the deformed samarium isotopes, we observe that as the excitation energy continues to 
increase K increases and reaches a plateau because of the contribution from rotational 
collective states before it decays in the vicinity of the shape phase transition energy. This 
behavior indicates that the rotational enhancement factor is correlated with the shape 
transition. This is reasonable since a spherical nucleus can no longer support rotational 
bands. 



CONCLUSION 



The microscopic calculations of statistical properties of nuclei at finite excitation ener- 
gies in the framework of the configuration- interaction shell model approach require very 
large model spaces. Such calculations have become possible using the SMMC method. 
We have extended SMMC to heavy nuclei and identified thermal signatures of the tran- 
sition from vibrational to rotational collectivity in families of rare-earth isotopes. 

We have found that the damping of vibrational and rotational collectivity as a function 
of excitation energy correlates with the pairing and shape phase transitions, respectively. 
Thus, the collective enhancement factor of the level density exhibits signatures of these 
phase transitions as a function of excitation energy. 
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